Large- iV /3-function for superconducting films in a magnetic field 
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Abstract By considering the large-N Ginzburg-Landau 
model, compactified in one of the spatial dimensions, we de- 
termine the beta-function and find an infrared stable fixed 
point for a superconducting film in presence of an external 
magnetic field, for dimensions 4 < D < 6. We find that 
this fixed point is independent of the film thickness (L). For 
D < 4, severe infrared divergences appearing under the form 
of divergence of the integral in the quantity a(D) defined af- 
ter Eq.(14), forbid an analogous analysis. For the physical 
dimension D = 3 we conclude that if a transition exists, it 
should not be a second order one. 
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It is usually assumed that it is a good approximation to 
neglect magnetic thermal fluctuations in the Ginzburg- 
Landau (GL) model, when applied to study high tem- 
perature superconductors. This is due to the fact that 
these materials have a very large Ginzburg parameter, 
typically k ~ 100. This problem has been investigated 
by a number of authors, both in its single component 
and in its N- component versions. A large account on 
the state of the subject and related topics can be found 
for instance in refs. [1-7] and other references therein. In 
particular in ref. [6] a large- N theory of a second order 
transition for arbitrary dimension D is presented and the 
fixed point effective free energy describing the transition 
is found. This author claims that it is possible that in the 
physical situation of N = 1, a mechanism of reduction 
of the lower critical domension could allow a continuous 
transition in D = 3. Here we investigate a confined ver- 
sion of the model studied in ref. [6]. We investigate the 
large- N limit of the Ginzburg-Landau model, the system 
being confined between two parallel planes a distance L 
apart from one another. Studies on confined field theory 
have been done in the literature since a long time ago. 
In particular, an analysis of the renormalization group in 
finite size geometries can be found in ref. [8] . This study 
is performed using a modified Matsubara formalism to 
take into account boundary effects on scaling laws. 

In this note we use an extended compactification for- 
malism in the framework of the effective potential, intro- 
duced in a recent publication [9], where the Euclidean 
massive (X(p 4 )D vector A-component model has been 
considered, and a study of the system subject to con- 
finement on a (i-dimensional subspace, d < D, has 
been performed. This allowed to generalize to a d- 
dimensional subspace, previous results in the effective 
potential framework for finite temperature and spacial 
boundaries. In particular it has been shown how a com- 



pactification mechanism applies to the study of spatially 
confined systems, for finite or zero-temperature field the- 
ories. This mechanism generalizes and unifies results 
from a recent work on the behaviour of field theories in 
presence of spatial planar boundaries [10,11], and previ- 
ous results in the literature for finite temperature field 
theory, as for instance in [12]. In the present work our 
aim is to show how the above mentioned formalism can 
be used to study in a field theoretical framework, effects 
associated to spacial confinement, for a model that allows 
a non-perturbative approach and is supposed to have a 
direct physical interpretation. We consider the vector 
A-component Ginzburg-Landau model in presence of an 
external magnetic field at leading order in the sys- 
tem being submited to the constraint of being confined 
between two parallel planes a distance L apart from one 
another. For k ^> 1 the Hamiltonian density of the GL 
model in an external magnetic field can be written in the 
form, 



H = \(V - ieA)<f>\ 2 + m 2 \<t>\' 



(i) 



where V x A = H and m 2 = a(T - T c ) with a > 0. 
This model Hamiltonian describes superconductors in the 
extreme type II limit. In the following we assume that 
the external magnetic field is parallel to the z axis and 
that the gauge A = (0, xH, 0) has been chosen. We will 
consider the model (1) with A complex components and 
take the large A limit at Nu fixed. If we consider the 
system in unlimited space, the field <f> should be written 
in terms of the well known Landau level basis, 
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where xi,p ,Pz( r ) are the Landau level eingenfunctions 
given by 
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with energy eigenvalues £/(|p|) = |p| 2 + (21 + \)lu + m 2 
and where Hi are the Hermite polynomials. In the above 
equations p and u are (D — 2)-dimensional vectors. 

Now, let us consider the system confined between 
two parallel planes, normal to the z-axis, a distance L 



apart from one another and use Cartesian coordinates 
r = (z, z), where z is a (D — 3)-dimcnsional vector, with 
corresponding momenta k = (fc z ,q), q being a (D — 3) 
-dimensional vector in momenta space. In this case, the 
model is supposed to describe a superconducting mate- 
rial in the form of a film. Under these conditions the 
generating functional of the correlation functions is writ- 
ten as, 
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(4) 



with the field <p(z,z) satisfying the condition of confine- 
ment along the z-axis, ip(z = 0, z) = ip(z = L, z). Then 
the field representation (2) should be modified and have 
a mixed series-integral Fourier expansion of the form, 

OC OO „ j „ 
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(5) 



where u> n = 2nn/L, the label / refers to the Landau levels, 
and the coefficients c„ and b(q) correspond respectively 
to the Fourier series representation over z and to the 
Fourier integral representation over the D— 3-dimensional 
z-space. The above conditions of confinement of the z- 
dependence of the field to a segment of length L, allow us 
to proceed with respect to the z-coordinate, in a manner 
analogous as it is done in the imaginary-time Matsubara 
formalism in field theory. The Feynman rules should be 
modified following the prescription, 
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(6) 



We emphasize that here we are considering an Euclidean 
field theory in D purely spatial dimensions, we are not 
working in the framework of finite temperature field the- 
ory. Temperature is introduced in the mass term of 
the Hamiltonian by means of the usual Ginzburg-Landau 
prescription. In the following we consider only the lowest 
Landau level I = 0. Thus, limiting ourselves to an analy- 
sis at the lowest Landau level (LLL), we obtain that the 
effective \<f>\ A interaction is given in momentum space and 
at the critical point by 



U(p) = 



1 + Nuoje-^^+Pl^ip)'' 



(7) 



where the single one-loop buble is given by, 
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The sum over n and the integral over k can be treated 
using the formalism developed in [9] . The starting point 
is an expression of the form, 
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(9) 



which, using a well-known dimensional regularization for- 
mula [8], can be written in the form, 



U = f(D, S )g s zf( S -^-l;a), 



(10) 



where f(D,s) is a function proportional to T(s — -j) 
and zf(s — -j;a) is one of the Epstein-Hurwitz zeta 
-functions defined by, 



Z c K (u;{ ai }) = {am\ + ... + a K n 2 K + c 2 )- u , 



ni,...,riK=—oo 



(11) 

valid for Re(u) > K/2 (in our case Re(s) > D/2). 
The Epstein-Hurwitz zeia-function can be extended to 
the whole complex s-plane and we obtain, after some 
manipulations [9], 



U = h{D,s) 
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(13) 



and K v are the Bessel functions of the third kind. 
Applying formula (12) to Eq.(8) the result is, 
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(14) 



where a(D) = J Q dx[x(l — x)]~ ^ . 

If an infrared stable fixed point exists, it would be pos- 
sible to determine it by an study of the infrared behaviour 



of the beta-function, i.e, in the neighbourhood of |p| = 0. 
Therefore we should investigate the above equations for 
|p| w 0. We start from (14), using an asymptotic formula 
for small values of the argument of Bessel functions, 

K v (z)K±r(u)(Zy (z ^ ; Re(u)>0), (15) 

which allows after some straighforward manipulations, to 
write Eq.(14) for |p| w in the form, 



S(p)«(27r) — 



2^a(i))r(^] |p|^- 6 + 



+L e-D 2 ^ T (D-6 H{D _ 6) 



(16) 

valid for odd dimensions D > 7, due to the poles of the F 
and C functions. We can obtain an expression for smaller 
values of D performing an analytic continuation of the 
Riemann zeta-function ((D—7) by means of the reflexion 
property of zeia-functions, 



C(*) = 



which gives, 
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(18) 



Inserting Eq.(18) in Eq.(16), and remembering the def- 
inition of a(D) after Eq.(14), we obtain an expression 
valid for 4 < D < 6, 



X(p)*A(D)\p\ D -s+B(D,L) 



(19) 



where, 



B(D,L) = n^L 6 - D T 



7-D 
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and 



A(D) = (2n)^2^T ) a(D). (21) 
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from which we get the infrared stable fixed point, 
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(25) 



Inserting (20) and (21) in Eq.(7) we have, 



Previous renormalization group calculations for mate- 
rials in bulk form in ref. [4] indicated a first order tran- 
sition for 4 < D < 6. This result has been obtained 
using an e = 6 — D expansion. The same conclusion 
of a first order transition is obtained in Ref. [1] with a 
large-N calculation. A large-N analysis and a functional 
renormalization group study performed in Refs. [5,6,13] 
leads to an opposite result to that of Ref. [1] , concluding 
for a second order transition in dimensions 4 < D < 6. 
The same conclusion is obtained in Ref. [7] . The authors 
in Ref. [5] claim moreover that the inclusion of fluctua- 
tions do not alter significantly the main characteristic of 
the system, that is, the existence of a continuous transi- 
tion into a spatially homogeneous condensate. Therefore 
for 4 < D < 6 our result for a film-like material (finite 
L) is the existence of an infrared stable fixed point, in 
agreement with those obtained for the material in bulk 
form (L = oo) in Refs. [5-7]. Moreover the fixed point is 
independent of the film thickness L. This could be inter- 
preted as pointing to a second order transition. However 
this should be taken as an indication not as a demon- 
stration of the existence of a continuous transition. As 
already discussed in [5,13], even if infrared fixed points 
exist none of them can be completely atractive, due to 
the large space of couplings. In this case the existence 
of an infrared fixed point as found in this paper, docs 
not assure the existence of a second order transition. For 
D < 4, severe infrared divergences appear under the form 
of divergence of the integral in the quantity a(D) defined 
after Eq.(14). We conclude that, for materials under the 
form of films, as is also the case for materials in bulk 
form, if there exists a phase transition for D < 4, in par- 
ticular in D = 3, it should not be a second order one. 
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1 + Nuoj L4(L>)| P p- 6 + B(D, L)) 
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Let us take as a running scale |p|, and define the dimcn- 
sionless coupling 



.9 = U a (p x =0,p y = 0, p)w|p 



D-6 



(23) 



where we remember that p is a D — 3-dimcnsional vector. 
Then we obtain straightforwardly the &eict-function for 
|p|»0, 
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